Phenomenological Analysis of pp and pp Elastic Scattering Based on Theoretical 

Bounds in High-Energy Physics 

S. D. CampofE 

Universidade Federal de Sao Carlos, campus de Sorocaba, 18052-780, Sorocaba, SP, Brazil. 

V. A. Okoroko^Q 

National Research Nuclear University "MEPhI", Kashirskoe Shosse 31, 115409 Moscow, Russia 

Considering the Froissart-Martin bound, Jin-Martin-Cornille bound and the optical theorem, we 
propose a novel parametrization for the total cross-section of proton-proton and antiproton-proton 
elastic scattering data. Using derivative dispersion relations we obtain the real part of the elastic 
scattering amplitude and thus the p parameter. Simultaneous fits to otot an d p are performed 
allowing very good statistical descriptions of the available data. Furthermore, predictions to 
and p at energies not used in the fit procedures are presented. For we obtain predictions at 
PsJ \ RHIC, LHC and future HC energies. 
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I. INTRODUCTION 



Nowadays, the study of high energy scattering is one of the most exciting topics in physics, and we have, as far 
we know, the exact theory of strong interactions, Quantum Chromodynamics (QCD), which can describe, based on 
Oh| the perturbative calculation scheme, the hadron-hadron interactions at short distances. However, interactions at 
large distances (near forward scattering) cannot be calculated even within perturbative approaches. Proton-proton 
(pp) and antiproton-proton (pp) elastic scattering are the simplest processes not yet explained in terms of a pure 
QCD description. The absence of such a description allows empirical analyses based on models (l|-[T3| or almost- 
independent models (some experimental data set possess theoretical bias) [15| whose main goal is to fit the available 
experimental data extracting useful information, contributing to the development of novel computational schemes. 
Several phenomenological models for high energy scattering are based on general principles of the Axiomatic Quantum 
Field Theory (AQFT), such as unitarity and crossing-symmetry, and on derived results as analyticity, and have proven 
to be successful in understanding or predicting the behavior of the hadronic scattering amplitude. 

Empirical parametrization for the total cross-section, differential cross-section, slope of the differential cross-section 
and others physical quantities have widely been used as a source of phenomenological modcl-indcpcndcnt indications of 
' the experimental data behavior. These approaches must be seen as a first step or attempt toward a formally rigorous 
phenomenological model-independent description of high energy elastic hadron scattering, embodying a predictive 
character. Moreover, important high energy theorems and bounds have been demonstrated in the last decades, 
providing rigorous formal constraints in the region of asymptotic energies, which cannot be disregarded in any reliable 
formalism, mainly related with phenomenological model-independent approaches. Of course, we must keep in mind 
that AQFT is a mathematical framework for the treatment and interpretation of relativistic quantum field theories 
and, in general, not all formal results based on AQFT results are satisfied by QCD. 

In this work, we present a novel parametrization based on formal high energy bounds taken from AQFT for the 
total cross-section of the elastic scattering amplitude. Specifically, we introduce a novel analytical parametrization for 
the total cross-section, inferred assuming some rigorous high energy bounds. The real part of the scattering amplitude 
is analytically evaluated through derivative dispersion relation. The crossing symmetry relates particle-particle and 
particle-antiparticle amplitudes. Fits to pp and pp total cross-section, p parameter are performed and may lead to a 
predictive and phenomenological model-independent approach, able to describe quite well the bulk of the experimental 
data. Furthermore, we present some predictions for the total cross-section for the LHC run at s/s = 7.0, 10, and 14 
TeV, for the future HC at y/s = 50, 100, 200 and 500 TeV and for the RHIC experiment at yfs = 200 and 500 GeV. 
It is important to stress that the parametrization proposed here allows us to obtain, under some physical conditions, 
the whole scattering function (beyond the forward direction). We intend to show these analyses elsewhere. 
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The paper is organized as follows. In section [TT1 the Froissart-Martin bound [ljl [lj}, Jin-Martin-Cornille bound 
fl8L Il9j and the optical theorem are used to propose a novel parameterization for pp and pp total cross-section and p 
parameter using only four fit parameters to each scattering reaction. In scction lHll wc explicitly obtain the real part of 
the forward elastic scattering amplitude and the p parameter. In section ITVl we show the individual fit results for <rtot 
and p. In section fVl we show the simultaneous fit results for <rtot & n d p- Section [VII presents a possible interpretation 
for the effective physical mechanism of this model comparing it with some models present in the literature. In section 
I VIII we present our Final Remarks with some criticism. 



II. TOTAL CROSS-SECTION 



The total cross-section is one of the most important physical observable in the elastic scattering. Its energy 
dependence has been a focus of intense theoretical interest since the establishment of QCD as the theory of strong 
interactions. However, at high energies there is no description based only on pure QCD formalism and the use of 
models with free fit parameters is the usual way to study pp and pp clastic scattering. In this section one obtains a 
novel paramctrization for the total cross-section based on well-established results taken from AQFT. 

We start considering the scattering amplitude -F(s, q 2 ) expressed as a function of two Mandelstam variables in the 
center-of-mass system, usually the energy squared s and the momentum transfer squared t = —q 2 . We can write 
F (s, q 2 ) in terms of its real and imaginary parts as 

F(s, q 2 ) = RcF(s, q 2 ) + ilmF(s, g 2 ), 

where the imaginary part of F represents the absorption in the scattering process and usually is related to the real 
part by dispersion relation (integral or derivative) that analytically connects pp and pp elastic scattering through 
crossing property. Therefore, if we know lm.F(s, q 2 ), then we are able to compute analytically the whole scattering 
amplitude (at q 2 = 0). 

The main goal of this section is to propose pp and pp total cross-section considering high energy bounds from formal 
results of the AQFT. A fundamental result derived from unitarity and analyticity asserts that high energy total cross- 
sections for the hadronic forward elastic scattering (at q 2 = 0) should be bounded according to Froissart-Martin 
bound [H [13 

7T 9 

ctot < — 9 In (s/s ), (1) 
mi 

where (Jtot = crtot (s) is the total cross-section of pp and pp and so is some energy initial value. Henceforth, we shall 
adopt So = 1 GcV 2 . This bound refers to an energy dependence of the total cross-section rising no more rapidly than 
In 2 s. As a remark, unfortunately, at present day we cannot be unambiguously discriminate between asymptotic fits 
of In s and In 2 s using high energy data (2p| - [2^ | . 

On the other hand, Cornille [19j improving the result of Jin and Martin fl8j proved the asymptotic lower bound 

atot > -jr. (2) 
Using both formal results ([1) and ([2]) we propose the following parametrization for the total cross-section 



<o P t = i + ^-r- 2 s , 



s 6 +-y y/2s a m% 

6 7T/3 



In 2 s. (4) 



tot s e+7 ^ V2 s a TO 2 



The first term in the right hand side of (J3j> and ((4]) will be responsible for the description of the data at low energies. 
Therefore, S and <5, 7 and 7 will be slightly different and at very high energies both terms must vanish. Yet, a and a 
also has a role at low energies allowing a more accurate description of the low energy data. Nevertheless, the second 
term in the right hand side of both parameterizations will be the leading one at high energies. Using the Pomeranchuk 
theorem, /3 and /3, a and a must tend to the same numerical value. 

As seen above, in the fitting process for the total cross-section data we expect a — > as well as 7 — > —6, and 
therefore, s 6+7 — ► 1 and s a — > 1 indicating the saturation of the Froissart-Martin bound. However, at the present 
day energy scale, we do not expect a = and so even really small values of a will give a slow increase with s for 
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the total cross-section when compared with the "full" Froissart- Martin bound. One consider here the a-parameter 
as some kind of "correction term" in the Froissart-Martin bound at low energies and it may represents some physical 
mechanisms that avoid the total cross-section saturation. 

As mentioned above, a may be treated as a correction in the Froissart-Martin bound at low energy. However, at 
(finite) high energies all the second term in the right hand side of Q and (J4| will behave as some power of s 



ir(3 



V2s a ml 



In 2 s 



(5) 



where c is a convenient constant and p the effective index. 

In order to apply the results obtained above to analyze the total cross-section data, we notice the important result 
assuring that if the Froissart-Martin is reached then the difference between pp and pp total cross-section goes as 

MM 
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which means that the difference can increases at most as In s and even in this case, 

jpp 



a 
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tot 

pp 

tot 



1, 



as s — > oo. This is the revised version of the Pomeranchuk theorem and in order to obey this formal result, we must 
impose the following constraints 



5 = 6, 7 = 7, a = a, /3 = /3. 



(0) 



Nevertheless, the energies achieved in the collision processes allow to release these constraints and therefore, as a 
first step, we adopt in the fitting procedure four independent parameters to pp and four to pp. 

III. REAL PART OF THE ELASTIC SCATTERING AMPLITUDE AND p PARAMETER 



In this section, we shall obtain the real part of the forward elastic scattering amplitude. In particular, we obtain 
p(s), the ratio of the real to imaginary part of the forward scattering amplitude in the high energy domain. 

The optical theorem connects the total cross-section to imaginary part of the elastic scattering at q 2 =0. On 
the other hand, since the imaginary part of the elastic scattering amplitude may be connected to the real part 
using dispersion relation one concludes that forward elastic scattering amplitude possess the same free parameters, 
except for the free-parameters coming from the subtraction terms of the dispersion relation, of the total cross-section. 
Connections between real and imaginary parts of the forward scattering amplitude have been widely investigated 
by dispersion relation in both integral and derivative forms (we use here the derivative form). Several papers have 
been devoted to calculations of the real parts of the pp and pp forward scattering amplitude, using a wide variety of 
dispersion relations and different representations for the energy dependence of the imaginary part of the scattering 
amplitude at q 2 = 0. Here, to obtain the real part of the elastic scattering amplitude, we use the approach adopted 
in [ill Hfil • We define two auxiliary functions (crossing between pp and pp) 



Im/± 1 



and relate these functions with the elastic scattering amplitude F pp /p p = / + ± /_ , where /+ = Re/+ + ilm/+ and 
/_ = Re/_ + iIm/_ and f± arc convenient amplitudes. In the forward direction, the derivative dispersion relation 
for even (+) and odd (-) amplitudes are expressed in terms of a tangent operator and in the case of one subtraction 
(equal to two subtractions in the even case) they are given by [27j 



Rej+M) = k 
s s 



tan 



2 dins 



Im/+( S ,0) 



Re/_(s,0) 



= tan 
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On the other hand, Bronzan, Kane and Sukhatme [28[ obtained derivative dispersion relations (without subtraction 
constant) using an integration parameter v. This parameter was used as free fit parameter or considered as constant 
over years p7j . 
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TT 
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tan 
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2 V dins 



lm/_(s,0) 



and Kang and Nicolescu [29| obtained derivative dispersion relations without v parameter (without subtraction con- 
stant) and 
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Avila and Menon [23] obtained the same result as Bronzan, Kane and Sukhatme without v parameter, which is 
equivalent to consider v = 1 in Bronzan, Kane and Sukhatme results. The odd case of Kang and Nicolescu can 



be obtained from Avila and Menon expressions replacing tan — ( 1 + — p — ) — » — cot — — p — and expanding the 

L 2 V dlns/J L 2 a In s J 

series around the origin in Avila and Menon derivative dispersion relations [271 ] . Furthermore, results of Avila and 
Menon are equivalent to Bronzan, Kane and Sukhatme considering the particular case v = 1. 



It is important to stress that replacement tan 



1 + 



din s 



cot 



2 dins 



may present some theoretical flow. 



As known, the first term of the cot(x) expansion presents the application of the inverse of the derivative operator. 
Mathematically speaking, let n a linear operator. Its inverse, n _1 , is by definition (H(x) is a non-zero analytic 
function) 

nTr 1 .^)] = u-^UHix)} = h(x). 

Furthermore, II -1 is also linear (30j . Of course, we consider x as an element of a linear space with a well-defined 
norm. In the context of the operatorial algebra, there exist II -1 if and only if to m > 



m\\H(x)\\ < \\UH(x) 



(7) 



To the especific case where II 



din / 



is always possible obtain m to satisfy the above result. However, we must 



take in to account the possible divergence result, i.e., 



din 



-H(s)=0, 



and therefore l/(d/dlns) — > oo. The above result may be viewed as a possible limitation to the functional form of 
the imaginary part of the forward amplitude: it cannot be constant in d/dlns. Of course, this result is valid only if 
the formal procedure presented above is correct. 

It has also been demonstrated by Fischer and Kolaf [3lM33j that at high energies the above even tangent operator 
can be replaced by its first order expansion. The odd representation also can be replaced by its first term [33j by 
using the notation of Kang and Nicolescu. However, Avila and Menon and Kang and Nicolescu representations are 
equivalent by the transformation showed above. Therefore, we retain here, for the sake of simplicity, only the first 
term of tangent series expansion in the Avila and Menon representation for both even and odd derivative dispersion 
relations. 

Re/+(s,0) _ k tt d Im/+(s,0) 
s s 2 d In s s 



5 



ggZ-M) = £ / 1 + _^_\ Im/_(a,0) 
s 2 \ dins J s 

This approach was recently used presenting good statistical results [HI . Of course, other terms of the tangent series 
expansion could be considered as perturbations of the first (leading) term. 
The real part of the scattering amplitude is given by 

ReF pp/pp (s,0) = Rc/+(a,0) ± Rc/_(s,0) 
s s s 

The most general versions of local quantum field theory lead to proving dispersion relation and, more generally, 
analyticity properties in two variables in a rather large domain, if one makes use of the positivity properties of 
the absorptive part of the scattering amplitude. Empirically, over many years, dispersion relation had always been 
consistent with measured data for energies reached by fixed target machines (e.g. pion nucleon scattering) or colliders 
(ISR and SppS colliders). Therefore, considering the above formalism we obtain (pp and pp) 



ReF pp (s,0) 



IT 

2 



(5(6 + 7) ir/3y/2 



s 6+ 7 
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ReF Pp ( Sl 0) 



6(6 + 7) 7T/V2 , 
7—^ 5—^- his 



a In 2 s 



'tot 



'tot J' 



(9) 



where the crossing property allows the simultaneous analysis to pp and pp scattering amplitude. Notice that pp and 
pp forward scattering amplitudes interact with each to other by means the crossing property. 

At high energies is expected that the imaginary part of the elastic scattering amplitude dominate the real part, 
i.e., the absorptive part will represents the whole scattering amplitude. Therefore, a simple measure of how fast the 
imaginary part increases is given by p, the ratio of the real to imaginary part of the forward scattering amplitude in 
the high energy domain 



and using ©, Q, 



RcF{s,0) 



ImF(s, 0) ' 

and @ we can explicitly write p(s) as (pp and pp) 



ffP(s) 
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In 2 , 



tot J 



(10) 



7T r pp _ pp 
^ l a tot 



tot J 



ln 2 s 



(11) 



The only free parameter in the above results is the subtraction constant k. It is not difficult to see that p(s) = 
p — > hiT 1 s as s — > co. In the next sections, we shall show both individual and simultaneous fit results to atot and p 
parameter. 



IV. INDIVIDUAL FIT RESULTS FOR a tot AND p 



Available experimental data on atot for proton-proton and antiproton-proton scattering were fitted by functions 
and (0J respectively. The database compiled by Particle Data Group has become a standard references and the 
corresponding computer readable files are used [34| . The statistic and systematic errors are added linearly in our 
analysis just as well as for some early analysis (see, for example (3514371 ]). 
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As seen the suggested parameterizations are valid only for s > s m - m , where s m [ n is some empirical low boundary (in 
GeV). During the analysis procedure we decrease the s mm value as much as possible in order to describe the wider 
energy domain with statistically reasonable fit quality. This logic is used in all our fit procedure below. Thus the fits 
have been made only at s > s m ; n for different values of low boundary. The fitting parameter values are presented in 
Table U for pp and in Table [TT1 for pp. Figure Q] shows the experimental energy dependence for fit for pp (Fig. [TJl) , pp 
(Fig. [TJd) and corresponding fits at ^/s m - m = 5 (red lines), ^/s m i n = 30 (blue lines). 

Paramctrization Q agree to experimental data very well for any ^/s m ; n under study. One needs to emphasize 
that ^ allows us to describe experimental points at ^/s m ; n = 3 very well. We do not see clear energy dependence 
of a and /3 within error bars but these parameters are very small and agree to expected values (Table A} . One can 
see the statistically acceptable fit qualities are obtained for y 7 ^min = 5 - f5 only for pp (Table P). Fit function (g]) 
agree to experimental points at qualitative level well for other low boundary values under study. But fit quality is 
statistically unacceptable at lowest ^/s m i n and some worse for y/s m [ Q = 20 — 30 1 . The a as well as /3 some decrease 
up to ^/s mm = fO and remain almost constants within errors at larger ^/s m ; n . On the other hand it seems small 
values of a, a agree with predicted behavior cr^^ oc In 2 s at high energies. One needs to emphasize the constraints 
([6]) are valid for fits at ^/s m j n > 15 within errors. 

In accordance with ([3]) and total cross-section energy dependence is described by the same functional form 
for pp and pp reactions. Thus one can make the quantitative analysis of common data ensemble for a^ t and (re- 
taking into account above results. The fitting parameter values are presented in Table UTTl for NN scattering. As seen 
the fit quality is statistically unacceptable for ^/s m ; n = 15 despite of equality of fit parameters for pp and pp. Our 
quantitative analysis demonstrates that experimental data for pp and pp total cross-section can be described by one 
function at ^/s m - m > 25 with statistically reasonable quality. Fig. [2] shows the experimental energy dependence for 
°tot an d corresponding fits for two various ^/s m ; n : fit at low energy limit for simultaneous study of proton-proton 
and antiproton-proton scattering ^s m - m = 15 is shown by red line, blue line corresponds to the fit with best quality 
a t \J s min = 25. 

If we fix the set parameters we are capable to obtain theoretical predictions for dtot at the RHIC, LHC and at the 
future HC energies. There are several predictions for the total cross-section covering a wide range, from 80 mb up to 
230 mb and at Tevatron (^s — 1.8 TcV) crtot is ill-determined, ctot = (71.42 ± 2.41) mb for E710 Collaboration and 
ctot = (80.03 ± 2.24) mb according to CDF Collaboration. Using y/s m i n = 5 as the minimum low boundary value 

.pp _ 

TeV. 

On the other hand, for y/s = 546 GeV (closest to the maximum RHIC energy), CDF Collaboration obtains 
°tot = (61.3 ± 1.0) mb, UA4 Collaboration otot = (61.9 ± 1.5) mb. Using ^Smin = 5 one obtains er^. = (61.93±6.17) 
mb for Vs = 546 GeV. 

At _Vi = 14 TeV, Block et al predict er^ t = (108.0 ± 3.4) mb [H] and Bourrely, Soffer and Wu predict of£ t = 103.5 
mb |39|. In Table HVl our predictions are shown for the pp total cross-section at LHC energy ranges, the future HC 
and at RHIC. To summarize, our predictions indicate cr to t — 100.44 mb at y 7 ^ = 14 TeV, crtot — 95.99 mb at \fs = 10 
TcV, and crtot — 91.35 mb at yfs = 7 TeV (possible LHC new energies). 

The extraction of the pp total cross-section from cosmic-ray experiments is based on the determination of proton- 
air production cross-section. The procedure is model-dependent [ITJ]. The Fly's Eye Collaboration 41 1 using the 



with statistically acceptable fit quality for pp data one obtains aPL = (73.22 ± 8.30) mb as a prediction for ^fs = 1 



Geometrical Scaling model obtain to pp at y 7 ! = 30 TcV crtot = (122 ± 11) mb. At same energy one obtains 
^tot = (110.69 ± 19.32) mb at y/s m - m = 3 as well as lowest energy boundary with statistically acceptable fit quality 
for pp scattering. On the other hand, Gaisser, Sukhatme and Yodh [42[ using the Fly's Eye result and the Chou-Yang 
prescription for the slope [43| obtained at ^fs = 40 TeV for pp interaction crtot = 175^? m b. At same energy one 
obtains crtot = (114.60 ± 20.55) mb. Of course, results from cosmic-rays presents a theoretical bias as noted by 
Nikolaev |44|. When the Akeno's result at y/s = 30 TeV (~ 92 mb) was corrected the corresponding crtot increase 
(amazing) on ~ 30 mb. 

Even at higher energies than that accessible in the present collider experiments, Baltrusaitis et al predicts a^ot ~ 
(120 ± 15) mb at ^fs — 30 TeV [il|. At same energy our prediction using the parameters obtained starting at 
= 3 is cr[£ t = (110.69 ± 19.32) mb. Finally, Block and Cahn [H] obtain at very high energies (y/s = 500 TeV) 
CT tot = (316.00 ± 5.46) mb and our prediction to the same energy is a^ t = (149.25 ± 32.85) mb, almost half. It is 
important to stress that the fit qualities for ([3]) and ((4]) are better for standard PDG paramctrization at y/s > 5 
GeV and suggested parameterizations allows to represent nucleon-nucleon crtot energy dependencies in wide energy 
domain. 



1 Values of fit parameters are equals at ^/s mul = 25 and for ^/s mul = 30 because of absence of experimental pp points in the range of 
collision energies y/s = 25 — 30 GeV. 
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The experimental data for ratio of the real to imaginary part of the forward scattering amplitude are some poorer 
than that for total cross-sections especially in high energy domain. Below in this paragraph individual fits under 
study only for p parameter energy dependence in pp and pp scattering. 

The fitting functions JTU]) and (| 1 1 1) are used for p pp (s) and for p^ p (s) respectively. 

There are two possible approaches for individual fit of p parameter for each interaction under considered, namely, i) 
the only one free parameter k is used in fitting functions (fTU|) and (fTT|) for p pp (s) and for p pp {s) but values of all remain 
parameters are fixed by the individual fits of total cross-sections at corresponding ^/s m j n ; and ii) all parameters in 
(fTU)) and (fTTj) are free, results of individual fits for a PP t and a PP t are not taken into account at all. The fit qualities are 
reasonable for any ^/s mm both for pp (Table [Vj part i) and pp ( Table IVTl part i) if the corresponding parameter values 
are fixed based on the above results for total cross-sections and only k parameter is free. Moreover, x 2 /n.d.f. are 
statistically acceptable at ^/s m ; n > 5 for pp and at some values of ^/s mm for pp in the framework of i) fit approach. On 
the other hand the precision of k parameter is some poor especially for pp interaction. The fit with all free parameters 
was studied also. The results for fit approach ii) are shown in Table |V] (part ii) for pp scattering and in Table IVII 
(part ii) for pp reaction 1 . For approach ii) high energy boundaries for individual fits of p have to be decreased because 
of large number of free parameters and small data samples for high energy domain. In this case the fit qualities are 
statistically acceptable for all s mm under study both for pp and pp interactions. 

As seen from Table fVl and from Table IVTl precision of fit parameters is some poorer than that for total cross-section 
fits for both approaches i) and ii). Fig. [3] shows the experimental energy dependence of p parameter for pp (Fig. 
[HK,b), for pp (Fig. [3fc,d) elastic scattering and corresponding fits at various values of low boundary. Fits at ^/s m ; n = 5 
are shown by red lines at all figures, fits at low boundary ^/s mm = 30 are available for i) approach only and shown by 
blue lines at Fig. [5Ji,c. The same color lines correspond to the ^s m - m = 20 (15) for pp (pp) reaction for ii) approach 
at Fig. [5}3,d. In the framework of i) approach fit curves show a noticeable increasing /decreasing of p parameter value 
for pp I pp interactions respectively for ^/s m ; n = 3 — 10 at high energies. The changing of absolute value of p is weaker 
at increasing of ^/s m ; n (Fig. [3^,,c) . At higher values of low boundary fit curves for pp demonstrate some maximum 
value in the y/~s ~ 100 GeV region but curves for pp clastic scattering show the much more flat behavior than that 
for pp reaction. Value of the maximum for pp reaction increases with increasing of low boundary for fitted energy 
domain at ^/s m ; n = 15 — 25. As seen from Fig. \Sjp the situation with increasing of low boundary is opposite in the 
framework of ii) approach for pp. The fit curves at Fig. [3ji show increasing of p parameter for pp interaction at any 
^/s mm under considered for ii) approach. It seems the model curves obtained in the framework of i) approach with 
one k free parameter show more reasonable and expected behaviour in the physics sense at ^/s m ; n > 15 at least than 
that in the framework of ii) approach. One can see a very small values for j3 at two highest values of ^ys m - m and for f3 
parameters in the case of individual fits of ppvlvv energy dependence. In the framework of ii) approach constraints ([6]) 
are valid for all fit parameters with exception of (5, S) pair at highest ^/s m in = 15 at which comparison is available 
for parameter values for pp and pp reactions. But one needs to emphasize that qualities for individual fits are well 
enough for ^/s mm = 3 even and suggested formulas (|10[) , (|11[) allow us to get a good agreement between model curves 
and experimental data in wide energy domain ^fs > 3 GeV. 

V. SIMULTANEOUS FIT RESULTS FOR a tot AND p 

The simultaneous fits were considered for various ensembles for scattering parameters. Results of simultaneous fit 
for {p^, fP p } by (fT0|) and (fTTj) with all free parameters are shown at Fig. and Fig. |4Jd for pp and pp, respectively. 
Significant increasing of fitted points allows us to study the approximation of p-parameter energy dependence up to 
highest y/s m [ n — 30. Numerical values of fit parameters are in the Table I VIII As seen fit qualities arc statistical 
acceptable at ^s m [ n > 3. Parameter errors are the same order or some smaller than that for individual fits at 
corresponding ^/s m ; n . Constraints (|6|) arc valid at ^/s mm > 15 for all parameters within two standard deviations at 
least. Approximation curves show smooth behavior at all y/s m i n . Some decreasing of p pp is observed at high energies 
(Fig. |4ji). This changing of pP p at high energies is most noticeable at ^/s m ; n = 15. Approximation curves are very 
close for p pp for 0.03 > yfs < 1 TeV at all ^/s m ; n under considered. Behaviour of fit curves for p pp (s) depends on value 
of ^/s m i n . At ^/s m j n = 15 and 30 corresponding curves demonstrate almost constant value p pp ~ 0.12 in t/s > 100 
GeV energy domain. But model curves show decreasing at high energies y/s > 1 TeV at all another values of low 
boundary. 



Values of fit parameters are equals at ^/s mm = 5 and for ^/s mm = 10 because of absence of experimental pp points in the range of 
collision energies ^fs = 5 — 10 GeV. 
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Simultaneous fits for {c^f °tot> P PP } by ©> ® anc ^ ar e shown at Fig. [5] for ^/s m ; n = 5 (red lines) and for 
^/s m ; n = 30 (blue line). The values of fit parameters are shown in Table [VIlTJ Suggested parameterizations allow us to 
obtain statistically acceptable fit qualities at ^s m [ n > 5 but fit quality is quite reasonable even at lowest ^s min = 3. 
The precisions are reasonable for all parameters. Constrains ([5]) are valid at ^/s m ; n > 15 for all parameters (with 
exception of <5, S at two lowest ^/s m j n = 15, 20 under considered in this case) within two standard deviations at least. 
Total cross-section curves for pp show a very close behavior both for ^/s m ; n = 5 and s m - m = 30. Model curve for o PP ^ 
at y'Smin = 5 shows significantly slow increasing for high energy domain than that at ^/s rri j n = 30. The approximation 
curve for p pp at ^/s m ; n = 5 demonstrates the opposite behaviour at collision energy increasing in high energy domain 
(y's > 1 TeV) than that the curve at highest low boundary for fitted energy domain. 

Results for simultaneous fits of {<J PP t , a PP tl pP p } by (J3|), ^ and (fTTj) are shown at Fig. |5]and in Table HXl Conclusions 
based on the numerical parameter values are similar them which have described above for Table IVIIII One can see 
that fit parameters have a close values both in Table IVIIII and in Table IIXI at the same value of lower boundary for 
energy domain at ^/s m j n > 3 with exception of fc-parameter. The total cross-section for pp shows significantly slow 
increasing in high energy domain at ^s m - m = 5 than that at y/s m - m = 30 for the ensemble of scattering parameters 
under study (Fig. [BJd) . Behaviour of model curves for p pp at high energies depends on value of low boundary. The 
curves show noticeable decreasing at ^s m - m = 5—15 for energy range y/s > 1 TeV. But fit curves demonstrate almost 
constant value p pp ~ 0.10 — 0.15 at higher values of low boundary for energy domain indicated above. 

Simultaneous fits of full set of global scattering parameters {af p , a PP , 

ppp^ppp) by ©, Q, (HDJ and CCQ) are shown 
at Fig. [7] at ^/s m ; n = 5 (red lines) and at ^/s m ; n = 30 (blue lines) . Numerical values of fit parameters are shown 
in Table [Xj Fit quality is statistically acceptable for ^s m ; n > 5. Moreover fit quality is reasonable at lowest 
value of y/s m [ n under study. Thus fit functions ([3]), (|4]), (fTO]) and (fTTj) allow us to obtain a reasonable description 
of experimental energy dependencies of global scattering parameters in wide energy domain. On the other hand 
one needs to emphasize that noticeable increasing (decreasing) of p pp (p pp ) value is observed for high energies at 
^/s m ; n < 10 (Fig. [7J;,d) especially at lowest boundary value ^/,s m ; n = 3. At larger ^s m - m there is a maximum in p(s) 
dependence with following much slower decreasing of p for energy domain ^fs > 1 TeV both in proton-proton and 
antiproton-proton interactions. Precisions are reasonable for all fit parameters at any low boundary values. These 
precisions are the same order or some smaller for this case than that for individual fits of scattering parameters. 
Constraints ([6]) are valid for all fit parameters at ^/s m i n > 20 within errors with exception for 5 and 5. The last two 
parameters have similar values, which are approached to each other with boundary value increasing. For S and 5 
parameters constraints ^ are valid at ^/s m - m indicated above within two standard deviations. Parameters a, a, /?, (3 
are almost constants within errors for ^/s m ; n > 15. 

The expected asymptotic behaviors are obtained for approximation curves for all scattering parameters under study 
at highest value of low boundary at least (Fig. [7|) . Curves for both pp and pp cross-sections shows a close behavior 
at all ^/Smin especially for proton-proton interactions. Corresponding curves for a pp demonstrate some differences 
for energies higher than 1 TeV at low y/s m [ n < 5 only. Difference between curves for p pp (p pp ) at several ^/s m i n are 
larger, especially for -y/s > 1 TeV energy domain. 

Predictions arc obtained for <r PP t , <r PP t at various energies based on the fit results of full set of global scattering 
parameters from Table [X] Total cross-section values are shown in the Table IXII These results are in agreement 
with predictions above from Table IPVl Uncertainties for these predictions are the same order or smaller for the case 
under considered than results, which were obtained for individual fits (Tabic llV|) . The smallest uncertainties are at 
^/s m ; n = 15. Perhaps, the one of the reasons for such feature is the optimal value of low boundary y/s m - m which 
allows to exclude the influence of spread of experimental points at low energies on uncertainties of fit parameters as 
well as to remain amounts of experimental data large enough for reasonable precision of fit. As seen from the Tabic 
IXII results for various ^/s m ; n are in a good agreement within errors. 

Our results coincide well with values of a PP t measured by UA4 and CDF Collaborations at ^/s = 546 GeV at any 
^/Sjnju, in particular, at y/s m - m = 15 which correspond to the best precision for predictions as well as statistically 
acceptable fit quality (see Table |X| . The similar situation at higher energy yfs = 1.8 TeV taking into account the 
wider dispersion of experimental values obtained by E710 and CDF Collaborations. For y/s m - m under considered 
above our predictions are in good agreement with Fly's Eye Collaboration [4l| data at v/s = 30 TeV but results in 
Table IXII are in qualitative agreement only with the estimation a pp t at y/s = 40 TeV [42J with taking into account 
errors. Total cross-section value at very high energy -y/s = 500 TeV significantly smaller than that from [45| as well 
as corresponding value obtained for individual fit results. 

It is important to stress that simultaneous fits are used because at high energies pp and pp elastic scattering 
amplitudes are (practically) the same function. Then, pp and pp experimental data tend to the same value at high 
energies and we expect that the absence of pp experimental data at some high energy value may be filled by presence 
of pp experimental data at (almost) same energy and vice-versa. Therefore, we expect that pp and pp simultaneous 
fit allows more reliable numerical results to the fit parameters. As can be viewed in Tables IIVI and |X] the errors in 
the parameters obtained from simultaneous fits are the same order or some smaller than those from individual fits 
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allowing a more precise description of the fitted physical quantity. 

Based on the above study one can conclude that suggested approximations allow to obtain fits (both individual and 
simultaneous) of experimental data for global scattering parameters with good statistical qualities for wide energy 
domain yfs > 5 GeV. Moreover this energy domain can be expanded as i/i > 3 GeV in some cases. The quality of last 
simultaneous fit at ^/s^\^ = 10 is better than that for Kang-Nicolescu model fit as well as for Donnachie-Landshoff 
fit of pp and pp experimental data [35j . Fit qualities obtained in this work for wide energy domains are close to the 
best fit qualities observed previously for various process in the framework of model with different number of poles 

[II 11 El E3. 

VI. COMPARISONS 

In this section, we compare the possible physical description of our model with some models present in the literature. 
There are several models but for the sake of the simplicity, we restrict ourselves to three models only. The first 
one is the famous and successful Donnachie-Landshoff model [48j based on a very simple parametrization structure 
evolving the exchange of reggeons and pomerons trajectories. The second one is the Kharzeev-Levin model j49T - [5l| 
based on N=4 Super-symmetric Yang-Mills theory and considering the contribution of D-instantons for multiparticlc 
production. Finally, the Godbole-Grau-Pancheri-Srivastava model p52l - [54| based on minijet contributions to the total 
cross-section. 



A. Donnachie-Landshoff Model 

One of the most successful models is the Donnachie-Landshoff (DL) parametrization [48| 

o%(a) = Xs< + Y pp s~\ aZ(s) = Xs* + Y»s~\ (12) 

where X, Y pp , Y pp , e, and r\ are fit parameters. The first term on the right hand side of the above expressions 
(X = 21.7 and e « 0.08) is associated with the pomeron contribution and has a common value to both pp and pp 
cross-section because the pomeron carries the quantum number of the vacuum and cannot distinguish particles and 
antiparticlcs. The second term (Y pp « 56.1, Y pp w 98.4 and r\ ~ 0.45) is associated with the reggeon trajectory, and 
it may be different to particles and antiparticles at low s. This term may be identified as resulting from p, w, /2, a>2 
exchange, for example. At s — > oo, the difference between pp and pp total cross-section parametrization can increases 
at most as Ins. However, the difference may vanish if a PP t = a PP t and, in any case, the ratio between both total cross 
sections tends to 1 as stated by the Pomeranchuk theorem. 

The DL model possesses clearly two regimes. Roughly speaking, at low energy the second term at the right hand 
side of Eqs. (|12[) is the leading one (the reggeon exchange dominates) and at high energy regime, the first term will 
be the leading one (the pomeron exchange dominates). This picture is quite similar to our model where the first 
term of the right hand side of ((3j) and (jlj is the sub-leading one and may be associated to the exchange of a reggeon 
trajectory with intercept apj(0) ~ 0.95. 

On the other hand, using ([5]) we may connect (at high energies) the behavior of the leading term of our model 
with the exchange of a pomeron trajectory. Adopting c = 21.7 as obtained in the DL model we obtain p < 0.06 
and therefore, ap(0) < 1.06. This value is slightly above 1 and may leads to the violation of the Froissart-Martin 
bound, but it does not occur for momenta lower than Planck scale. In this way, the pomeron exchange mechanism 
may ensures the preservation of unitarity. 

However, comparisons done above are not exactly correct. In our model, at low energies, reggeon and pomeron 
trajectories interacts with each other mixing the intercept contributions. Therefore, the reggeon intercepts showed 
here present pomeron contributions. Indeed, f% trajectory, for example, has the quantum number of the vacuum 
as the pomeron. There may be contributions from pomerons at low energies, preventing the complete separation 
between pomerons and reggeons. We expect, at high energies, the divergence of reggeon and pomeron trajectories 
from each other. Then, the influence of reggeons will ceases and only the pomeron exchange will be the main physical 
mechanism restoring the validity of Pomeranchuk theorem. 

B. Kharzeev-Levin Model 

The Kharzeev-Levin (KL) model [49T - l5l| is based on the Super-symmetric Yang-Mills theory in N=4 (N=4 SYM). 
Considering contributions of D-instantons in AdS§ bulk space, they argue that D-instantons coupled to dilatons and 



10 



axions are responsible for multiparticlc production in strongly coupled N=4 SYM. As a consequence, the cross-section 
increases with the energy 

^inel oc s A ~ 2 , 

where A depends on the convenient choice of parameters [49j-i5~D] . This approach is interesting since it suggests that 
topological effects may be important in high energy collisions. Furthermore, they argue that the weakly coupled 
graviton generates the elastic amplitude and the correspondent part of the total cross-section. On the other hand, 
D-instanton induced interactions of dilatons and axions are responsible for the multiparticlc production processes. 

This formalism predicts that the cross-section will increase with the energy due to D-instanton mechanism. This 
mechanism is induced by interactions of dilatons and axions responsible for the multiparticlc production process. The 
pomeron intercepts predicted in KL model depends on the value of A (and (J\ ne \ <x s A ~ 2 ). They found A sa 4 and 
consequently, the slope a' P rts 0.5 GeV~ 2 . 

However, this value of A may imply in a strong unitarisation scheme to restore the Froissart-Martin bound. As 
shown above, the pomeron exchange may be the main mechanism responsible for restore the unitarisation. Indeed, 
they obtain a very high pomeron intercept, possibly indicating the exchange of a new soft pomeron trajectory family 
(A is not unique). 

The model proposed here possesses a very smooth increase (at high energies) when compared with the standard 
theoretical result atot ~ 60 In 2 s. The KL model may explains these differences in the behaviors by the mechanism of 
multiparticlc production: at high energies the multiparticlc mechanism acts, increasing the hadron cross-section and 
when the saturation limit approaches the pomeron exchange has become physically important restoring the unitarity 
scheme turning the increase of the total cross-section smooth as obtained. Moreover, we may suppose a mixed state 
" reggcon+pomeron" as being the responsible for the observed differences between pp and pp total cross-section. 



C. Godbole-Grau-Pancheri-Srivastava Model 



The study of QCD minijet contribution is the main goal of the Godbole-Grau-Pancheri-Srivastava (GGPS) model 
[52Tl5"H . This model includes a re-summation of soft gluon radiation which they argue it is necessary to tame the fast 
high energy rise typical of a purely perturbative minijet model. 

The minijet formalism predicts a very fast rise to <7j ct sa s 3 at high energies. To restore the total cross-section 
they suppose that soft gluon emissions from the colliding partons are the physical attenuation mechanism. To restore 
the Froissart-Martin bound, they obtain a "new" formulation of the Froissart-Martin bound 

o-tot ~ (eln S ) 1/n 

where e fixes the asymptotic rise and 1/2 < n < 1 modulates the infrared behavior of a s f52Tj54j ] . 

Pomeron exchange is believed to be the exchange of a system of gluons. Therefore (ultimately), the pomeron 
exchange will be the Froissart-Martin bound restoration mechanism as in the KL model. As seen above, the pomeron 
exchange has a very important role if we want to understand the taming mechanism behind no saturation effect in 
our model. 



VII. FINAL REMARKS 



The model proposed here to describe the hadron-hadron total cross-section and the p parameter is basically based 
on two rigorously bounds obtained in the framework of AQFT. Namely, the Froissart-Martin bound and Jin-Martin- 
Cornille bound. Roughly speaking, in this model the like- Froissart-Martin bound will control the high energy behavior 
of the total cross-section, and like- Jin-Martin- Cornille bound will defines its behavior at low energies. However, the 
analysis of our model is not so simple, since we suppose an attenuation mechanism in the Froissart-Martin bound 
~ s~ a In 2 s which mix these two bounds, especially at low energies. 

The statistical description of ctot an d p experimental data is acceptable. The quality of some simultaneous fits is 
better than that for Kang-Nicolescu model fit as well as for Donnachie-Landshoff fit of pp and pp experimental data. 
The qualities of the fits obtained in this work for wide energy domains are closer to the best fit qualities previously 
observed for several processes in the framework of the model with a different number of poles. Therefore, based on this 
statistical fact, it is possible perform predictions in the energy variable and the results can be compared with others 
present in some models and parameterizations in available literature. Results viewed here shown some agreement 
among themselves. 
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Comparisons with some models may allow some physical understanding of what is happening here. In the context of 
DL model, we may understand the the index in the low energy term of our model as representing a reggeon exchange 
trajectory, and the index in the high energy term - as a pomeron exchange trajectory. Yet, the picture viewed here 
seems to be not so simple. The mixed description proposed, especially at low energies, prevents the identification 
of each index contribution separately (or each trajectory separately) and, therefore, we cannot treat this low energy 
term as being a pure reggeon exchange. 

In the context of KL and GGPS models, the pomeron exchange is the responsible for the attenuation of the rise 
of the total cross-section. However, in both models the saturation of the total cross-section is quickly reached and 
therefore, it is necessary the exchange of not only one pomeron but a whole family. In our model, this attenuation 
mechanism may be explained by the interaction of pomerons and reggeons that tames the rise of the total cross- 
section and this mixed state of index (trajectories) may be responsible for the observed differences between pp and 
pp behaviors. 

On the other hand, the model presented here does not saturate the total cross-section. Then, we may suppose that 
this is due to the exchange of more than only one pomeron trajectory since the increase of experimental data set 
in the fit procedures may give rise to different values of the effective index p and each value may represent a novel 
pomeron trajectory. Therefore, novel experimental results to total cross section can be added to the original data set 
allowing more precise values to the effective index p through fit procedures. These trajectories may tend to a limiting 
value, i.e, pi — > pi where i represents a pomeron trajectory and I represents the limiting trajectory [ap i (0) — > op,(0)). 
If this pomeron trajectory limit exist it may represent a definite (or not, since a new physical attenuation mechanism 
may be discovered) bound to the rise of total cross section. This behavior at very high energies is one of the puzzles to 
be solved at the future HC because the applied mechanisms used to tame its very fast rise, we hope, will be verified. 
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3 
5 

10 
15 
20 
25 
30 



0.028±0.008 
0.020±0.010 
0.010±0.012 
0.010±0.008 
0.019±0.003 
0.031±0.002 
0.021±0.012 



P x 10 3 
3.9±0.3 
3.4±0.5 
2.9±0.6 
2.9±0.3 
3.6±0.1 
4.6±0.1 
3.7±0.9 



-7 

5.882±0.005 
5.890±0.010 
5.910±0.020 
5.906±0.005 
5.876±0.004 
5.821±0.009 
5.870±0.040 



51.8±0.5 
50.8±1.5 
48.5±2.8 
48.8±1.0 
54.8±1.1 
69.5±4.2 
55.7±10.1 



X 2 /n-d-f 
107.6/139 
74.9/107 
46.7/69 
34.2/58 
30.1/45 
15.9/33 
15.1/30 



TABLE I: Parameters for fitting by ([3]) of ofj?^ energy dependence. 



3 
5 

10 
15 
20 
25 



a p x 10 3 -7 5 x 2 /n-d-f 

0.085±0.008 9.7±0.3 5.664±0.006 130.9±1.8 160.6/69 
0.059±0.006 6.9±0.5 5.749±0.013 99.8±3.8 35.7/55 
0.040±0.013 5.0±1.0 5.808±0.034 80.8±9.0 22.4/25 
0.005±0.026 2.7±1.2 5.895±0.043 58.1±8.5 16.3/20 
0.013±0.032 3.1±1.8 5.875±0.067 63.6±14.8 15.2/13 
0.023±0.031 3.8±2.1 5.836±0.094 77.1±28.1 14.8/11 



TABLE II: Parameters for fitting by Q of cr?]L energy dependence. 



^(TeV) 

0.2 
0.5 
7 

10 

14 

30 

40 

50 
100 
200 
500 



51.95 ±4.32 
60.46 ±6.14 
91.4 ±13.7 
96.0 ±15.0 
100.4±16.3 
110.7±19.3 
114.6±20.6 
117.7±21.5 
127.2±24.7 
136.7±28.1 
149.3±32.9 



51.06±6.62 
59.66±9.15 

92.6 ±19.7 

97.7 ±21.5 
102.7±23.3 
114.4±27.7 
118.9±29.4 
122.4±30.8 
133.7±35.5 
145.3±40.4 
160.8±47.5 



10 

52.09±9.09 
61.0±12.2 
96.7 ±26.0 
102.6±28.5 
108.3±30.9 
122.0±37.0 
127.4±39.5 
131.6±41.5 
145.4±48.2 
159.8±55.4 
179.8±66.1 



15 

<t ( mb ) 
51.42 ±4.58 

60.29 ±6.55 

96.1 ±15.3 

102.0±16.8 

107.7±18.4 

121.4±22.3 

126.8±23.9 

131.1±25.2 

144.9±29.4 

159.3±34.1 

179.4±41.0 



20 

52.42±1.73 
61.82±2.41 
98.04±5.53 
103.75±6.09 
109.26±6.64 
122.20±8.01 
127.23±8.56 
131.18±9.00 
143.7±10.5 
156.5±12.1 
173.9±14.3 



25 

52.84±1.74 
62.64±2.06 
97.92 ±3.98 
103.13±4.33 
108.09±4.67 
119.42±5.50 
123.71±5.83 
127.04±6.09 
137.35±6.94 
147.56±7.84 
160.79±9.10 



30 

52.0±12.3 
61.3±15.5 
96.8 ±29.7 
102.3±32.2 
107.7±34.6 
120.2±40.5 
125.0±42.9 
128.8±44.7 
140.7±50.9 
152.9±57.4 
169.3±66.7 



0.546 
1.8 



62.10±2.94 
70.37±3.80 



61.93±6.17 
73.22±8.30 



61.9±14.0 
74.9±19.8 



"tot 



(mb) 
61.4±26.8 
76.2±38.8 



60.8±35.6 
75.1±50.3 



61.3±36.9 
75.6±50.7 



TABLE III: Predictions for at various energies based on the individual fits of energy dependence of proton-proton and 
proton-antiproton total cross-sections. 



\/ s min 


ajvjv 


Pnn x 10 3 


15 


0.042±0.008 


5.3±0.7 


20 


0.034±0.011 


4.6±0.9 


25 


0.011±0.017 


2.9±1.0 


30 


0.010±0.003 


2.8±0.1 



-'JNN 



Snn x /n-d.f 



5.804±0.028 73.7±7.7 519.3/82 
5.837±0.036 64.1±8.8 196.9/62 
5.912±0.040 47.3±7.0 91.1/48 
5.915±0.004 46.6±1.1 90.8/45 



TABLE IV: Parameters for fitting by ((3]) of energy dependence. 
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Parameter 








V s min 








3 


5 


10 


15 


20 


25 30 


i) 


k, mbarn 


39 ±6 


-0.5 ± 17 


-239 ± 36 


15 ±60 


-115 ± 104 


571 ± 239 648 ± 301 


X 2 /n.d.f. 


146/80 


93.9/63 


39.1/47 


32.3/33 


13.8/22 


1.53/9 0.16/6 


ii) 


k, mbarn 


32 ±8 


150 ± 19 


514 ± 2 


1202 ± 2 


40.0 ± 1.4 




Q 


0.54 ±0.05 


0.55 ± 0.09 


0.53 ± 0.06 


-0.20 ±0.11 


-0.15 ±0.12 




f3 x 10 3 


4.4 ±0.5 


6±5 


0.8 ±0.2 


0.006 ± 0.003 


0.0011 ± 0.0002 




-7 


5.31 ±0.05 


5.21 ± 0.10 


5.21 ±0.06 


4.4 ±0.1 


4.10 ±0.18 




S, mbarn 


40 ±20 


53 ±62 


50 ±27 


141 ± 55 


129 ± 12 




Q 


0.59 ±0.17 


0.7 ±0.3 


0.7 ±0.5 


-0.18 ±0.07 


-0.13 ±0.11 




(3 x 10 3 


0.66 ±0.07 


2±2 


0.057 ± 0.005 


0.01 ±0.01 


0.0022 ± 0.0003 






5.12 ±0.08 


5.04 ±0.07 


4.963 ± 0.003 


5.010 ±0.006 


5.041 ±0.013 




S, mbarn 


51 ±3 


181 ±38 


733 ±6 


1443 ± 22 


40 ±3 




X 2 /n.d.f. 


75.3/72 


63.9/55 


38.9/39 


26.0/25 


12.4/14 





TABLE V: Parameters for individual fitting by (|10[1 of p pp energy dependence. Part i) corresponds to approach with one free 
parameter k and part ii) shows results for approach with all free parameters in (|10[) . 



Parameter 








V s min 








3 


5 


10 


15 


20 25 


30 


i) 


k, mbarn 


109 ± 25 


65 ± 144 


38 ± 131 


-436 ± 172 


-87 ±416 1172 ± 1197 i 


S84± 1197 


X 2 /n.d.f. 


44.4/17 


20.8/11 


16.6/11 


2.73/10 


8.21/6 11.7/4 


7.60/4 


ii) 


k, mbarn 


-0.3 ±0.2 


15 ±3 




18 ±3 






a 


-0.4 ±0.2 


0.078 ±0.016 




-0.03 ±0.10 






/3 x 10 6 


0.13 ±0.10 


4.5 ±0.9 




3±3 






-7 


4.0 ± 1.3 


4.8 ±0.4 




4.55 ±0.17 






S, mbarn 


-5 ±4 


10 ± 11 




61.5 ±48.2 






a 


-0.4 ±0.2 


0.056 ±0.015 




-0.05 ±0.10 






[3 x 10 6 


0.06 ± 0.05 


3.4 ±0.6 




2.0 ± 1.3 






T7 


5.52 ±0.02 


5.07 ±0.06 




5.12 ±0.12 






S, mbarn 


0.5 ±0.5 


13 ±3 




13 ± 10 






X 2 /n.d.f. 


9.67/9 


1.84/3 




1.78/2 







TABLE VI: Parameters for individual fitting by (|lip of p pp energy dependence. Part i) corresponds to approach with one free 
parameter k and part ii) shows results for approach with all free parameters in (|lll) . 



Parameter 


















3 


5 


10 


15 


20 


25 


30 


k, mbarn 


11.7 ±3.5 


21 ±7 


-23 ± 19 


109 ± 42 


-242 ± 68 


-174 ± 147 


-200 ± 189 


a 


0.05 ±0.02 


0.051 ± 0.007 


0.038 ± 0.006 


0.06 ± 0.03 


0.045 ± 0.006 


0.027 ±0.008 


0.05 ±0.03 


/3 x 10 3 


0.92 ±0.18 


0.78 ±0.03 


1.12 ±0.06 


1.10 ±0.15 


1.06 ±0.07 


1.05 ±0.09 


1.3 ±0.7 


-7 


5.84 ±0.02 


5.844 ±0.016 


5.866 ± 0.009 


5.68 ±0.06 


5.876 ±0.013 


5.84 ±0.02 


5.80 ±0.13 


S, mbarn 


10.05 ±0.49 


8.2 ±0.8 


13.3 ±0.8 


22 ±3 


10.6 ± 1.0 


16 ±2 


18 ± 12 


a 


0.061 ±0.005 


0.065 ± 0.006 


0.043 ± 0.005 


0.05 ± 0.02 


0.045 ± 0.006 


0.023 ± 0.008 


0.03 ±0.02 


(3 x 10 3 


1.18 ±0.10 


1.05 ±0.05 


1.28 ±0.08 


0.91 ±0.09 


1.11 ±0.09 


0.94 ±0.11 


0.9 ±0.3 


r7 


5.73 ±0.02 


5.705 ± 0.004 


5.802 ± 0.009 


5.75 ±0.05 


5.855 ±0.013 


5.887 ±0.016 


5.90 ±0.10 


S, mbarn 


19.7 ±2.5 


18.7 ±0.6 


20.6 ± 1.1 


16.2 ± 1.9 


12.2 ± 1.1 


12.9 ± 1.5 


11 ±6 


X 2 /n.d.f. 


86.6/90 


66.3/67 


42.1/51 


29.8/36 


12.8/21 


1.18/6 


0.25/3 



TABLE VII: Parameters for simultaneous fitting by (|10[) and (|ll|l of p energy dependencies for proton-proton and antiproton- 
proton collisions respectively. 
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Parameter 








V s min 










3 


5 


10 


15 


20 


25 


30 


k, mbarn 


42.0 ±6.5 


-52 ± 18 


-131 ±42 


32 ±20 


107.6 ± 74.5 


671.4 ±426.5 


621 ± 454 


a 


0.016 ±0.007 


0.018 ±0.003 


0.011 ± 0.004 


0.017 ±0.005 


0.025 ±0.010 


0.024 ±0.015 


0.022 ± 0.008 


P x 10 3 


3.4 ±0.3 


3.44 ± 0.07 


3.0 ±0.1 


3.42 ±0.15 


4.1 ±0.6 


4.0 ± 1.3 


3.8 ±0.6 


-7 


5.888 ± 0.004 


5.889 ±0.001 


5.903 ± 0.004 


5.885 ± 0.003 


5.85 ±0.02 


5.85 ±0.07 


5.865 ± 0.031 


8, mbarn 


51.3 ±0.4 


51.2 ±0.4 


49.1 ± 0.9 


52.6 ± 1.1 


60 ±6 


61 ± 17 


57 ±8 


a 


0.084 ± 0.003 


0.056 ± 0.005 


0.042 ± 0.011 


0.013 ± 0.003 


0.026 ± 0.018 


0.023 ± 0.012 


0.023 ± 0.011 


P x 10 3 


9.6 ±0.3 


6.4 ±0.5 


5.15 ±0.90 


3.13 ±0.13 


4.0 ± 1.3 


4± 2 


4±2 


T7 


5.669 ± 0.005 


5.758 ±0.012 


5.80 ±0.03 


5.879 ± 0.002 


5.84 ± 0.05 


5.84 ±0.10 


5.84 ± 0.09 


S, mbarn 


129.3 ± 1.7 


97 ±3 


83 ±9 


61.8 ± 1.2 


74 ± 14 


77 ±34 


77 ±31 


X 2 /n.d.f. 


392/288 


182/225 


109/141 


81.1/111 


60.2/80 


31.7/53 


30.1/47 



TABLE VIII: Parameters for simultaneous fitting by ©, @ and ()10l) of G lnt an< ^ P VV ener gy dependencies. 



Parameter 








\/ s mm 










3 


5 


10 


15 


20 


25 


30 


k, mbarn 


102 ± 20 


21 ± 17 


-91 ± 39 


-399 ± 188 


-585 ± 208 


-230 ± 123 


-270.5 ± 109.7 


a 


0.035 ± 0.006 


0.026 ± 0.004 


0.011 ±0.003 


0.018 ± 0.004 


0.017 ±0.005 


0.018 ±0.004 


0.017 ±0.003 


P x 10 3 


4.2 ±0.3 


3.77 ±0.13 


2.95 ±0.07 


3.38 ±0.11 


3.40 ±0.17 


3.31 ±0.15 


3.35 ± 0.08 


-7 


5.878 ± 0.005 


5.884 ± 0.002 


5.906 ± 0.002 


5.899 ± 0.003 


5.884 ± 0.002 


5.880 ± 0.004 


5.889 ± 0.005 


5, mbarn 


52.15 ±0.44 


51.6 ±0.4 


48.6 ±0.7 


51.5 ± 1.1 


53.0 ± 1.3 


54.1 ± 1.2 


51.6 ± 1.8 


a 


0.082 ± 0.003 


0.058 ± 0.002 


0.036 ± 0.002 


0.014 ±0.003 


0.012 ± 0.003 


0.016 ±0.013 


0.013 ±0.005 


P x 10 3 


9.4 ±0.3 


6.63 ±0.12 


4.69 ±0.13 


3.15 ±0.11 


3.09 ±0.11 


3.3 ±0.7 


3.2 ±0.2 


r7 


5.669 ± 0.005 


5.750 ± 0.002 


5.818 ±0.003 


5.880 ± 0.003 


5.878 ± 0.006 


5.865 ±0.025 


5.870 ± 0.007 


S, mbarn 


129.7 ± 1.7 


99.6 ±0.8 


78.1 ± 1.3 


61.3 ± 1.4 


63 ±3 


67 ±7 


65 ±5 


X 2 /n.d.f. 


288/225 


115/173 


72.3/105 


53.4/88 


46.6/64 


31.7/48 


30.9/45 



TABLE IX: Parameters for simultaneous fitting by Q, (|4j) and (|11[) of u P ^ t , cr?? t and p pp energy dependencies. 



Parameter 


3 


5 


10 


V s min 
15 


20 


25 


30 


k, mbarn 


48 ±6 


-51 ± 18 


-121 ±42 


-64 ± 35 


32 ± 14 


470 ± 371 


400 ± 149 


a 


0.025 ± 0.006 


0.021 ± 0.008 


0.014 ±0.010 


0.008 ± 0.004 


0.022 ± 0.009 


0.016 ±0.010 


0.013 ±0.012 


P x 10 3 


3.7 ±0.3 


3.5 ±0.4 


3.1 ±0.5 


2.9 ±0.1 


3.8 ±0.6 


3.3 ±0.6 


3.1 ±0.7 


-7 


5.883 ± 0.004 


5.888 ± 0.009 


5.900 ± 0.014 


5.908 ± 0.002 


5.87 ±0.02 


5.90 ±0.03 


5.90 ± 0.03 


5, mbarn 


51.8 ±0.4 


51.3 ± 1.2 


50 ±2 


48.4 ±0.9 


56 ±5 


52 ±5 


49 ±5 


a 


0.080 ± 0.003 


0.053 ± 0.005 


0.038 ± 0.010 


0.014 ±0.003 


0.026 ±0.014 


0.022 ±0.016 


0.019 ±0.013 


P x 10 3 


9.1 ±0.3 


6.2 ±0.4 


4.9 ±0.7 


3.18 ±0.12 


4.0 ±0.9 


3.8 ±0.1 


3.6 ±0.8 


r7 


5.676 ± 0.005 


5.763 ±0.011 


5.81 ±0.02 


5.874 ± 0.002 


5.84 ± 0.03 


5.84 ± 0.04 


5.85 ± 0.03 


S, mbarn 


127.3 ± 1.6 


96 ±3 


81 ±7 


63.0 ± 1.2 


74 ± 10 


75 ± 14 


72 ±9 


X 2 /n.d.f. 


423/306 


186/237 


113/153 


87.3/122 


62.2/87 


33.1/58 


31.4/52 



TABLE X: Parameters for global simultaneous fitting by ©, ©, JT0]) and (fTTj) of a PP t , af P t , p vv and p pp energy dependencies. 



y^(TeV) 

0.2 
0.5 
7 

10 
14 

30 

40 

50 
100 
200 
500 



51.34 ±3.80 
59.87 ±5.30 
91.5 ±11.4 
96.3 ±12.4 
101.0±13.4 
111.8±15.8 
115.9±16.8 
119.1±17.6 
129.3±20.2 
139.6±22.9 
153.2±26.7 



51.53±5.33 
60.25±7.34 
93.3±15.8 
98.5±17.2 
103.5±18.6 
115.1±22.1 
119.6±23.5 
123.2±24.6 
134.4±28.3 
145.8±32.2 
161.1±37.9 



10 

51.55±7.12 
60.32±9.69 
95.0 ±20.8 
100.6±22.8 
106.1±24.7 
119.0±29.5 
124.0±31.5 
128.0±33.0 
140.8±38.2 
154.0±43.8 
172.3±52.0 



15 



'tot 



(mb) 
52.37 ±1.93 
61.68 ±2.84 
99.49 ±7.10 
105.73±7.89 
111.84±8.68 
126.5±10.7 
132.3±11.5 
136.9±12.1 
151.7±14.3 
167.2±16.8 
189.0±20.4 



20 

52.66±7.82 
62.1 ±10.3 
97.6 ±20.7 
103.2±22.5 
108.5±24.3 
120.9±28.6 
125.7±30.3 
129.4±31.6 
141.2±36.1 
152.3±40.9 
169.4±47.6 



25 

53.69±9.60 
62.7 ±12.0 
97.7 ±23.1 
103.3±25.0 
108.8±27.0 
121.6±31.9 
126.7±33.9 
130.6±35.4 
143.2±40.7 
156.3±46.3 
174.1±54.5 



30 

51.7 ±10.6 
60.6 ±13.7 
96.4 ±27.7 
102.2±30.2 
107.8±32.6 
121.2±38.7 
126.4±41.2 
130.6±43.2 
143.8±49.7 
157.7±56.8 
176.7±67.1 



0.546 
1.8 



64.29±3.02 
71.28±3.92 



62.44±5.24 
74.54±7.11 



62.4 ±10.7 
75.7 ±14.8 



<ot M>) 
61.17±2.61 
75.60±3.89 



62.1 ±15.4 62.2 ±11.9 62.2 ±14.8 
76.1 ±21.6 76.8 ±17.4 77.0 ±20.8 



TABLE XI: Predictions for cr^ot a t various energies based on the simultaneous fits of energy dependence of proton-proton and 
proton-antiproton total cross-sections. 
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FIG. 2: (color online) Total cross-section energy dependence for nucleon-nucleon collisions. The red line corresponds to the fit 
at y/^mhi = 15, the blue line - at ^/s m ; n = 25. Experimental points are indicated as • (■) for pp (pp). Experimental data are 
from 1341 1 . 
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FIG. 3: (color online) The p parameter energy dependence for proton-proton (a,b) and antiproton-proton (c,d) collisions. 
Results obtained in the framework of i) approach for individual fit are shown on (a,c) and for approach ii) - on (b,d). The 



red lines correspond to the individual fits at ^fs 



5 for any approaches. Blue lines correspond to the individual fits at 



^/s m - m — 30 for i) approach (a,c) and for approximation curves obtained in the framework of ii) approach at ^/s m ; n = 20 (15) 
for pp (pp) respectively (b,d). Experimental data are from 
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FIG. 4: (color online) The p parameter energy dependence for proton-proton (a) and antiproton-proton (b) collisions and results 
of simultaneous fits of these two parameters. The red line corresponds to the fit at ^/s m j n = 5, the blue line - at ^/s m ; n = 30. 
Experimental data are from [34| . 
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FIG. 5: (color online) The a"tot ener gy dependence for proton-proton (a), antiproton-proton (b) collisions, p pp parameter vs 
collision energy (c) and results of simultaneous fits of these parameters. The red line corresponds to the fit at ^/s m ; n = 5, the 
blue line - at ^/s m ; n = 30. Experimental data are from [3l[ . 
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FIG. 6: (color online) The a"tot ener gy dependence for proton-proton (a), antiproton-proton (b) collisions, p pp parameter vs 
collision energy (c) and results of simultaneous fits of these parameters. The red line corresponds to the fit at ^/s m ; n = 5, the 
blue line - at ^/s m ; n = 30. Experimental data are from [3l[ . 
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FIG. 7: (color online) The crtot ener gy dependence for proton-proton (a), antiproton-proton (b) collisions, p vv (c), p vv (d) 
parameters vs collision energy and results of simultaneous fits of all fo ur p arameters. The red line corresponds to the fit at 
^/s m ; n = 5, the blue line - at ^/s m ; n = 30. Experimental data are from [33 ]. 



